For the problem of low reservoir prediction accuracy near the fault, the fault area is divided into “stable area” and “complex area”. By precisely interpreting the “stable area” of faults and establishing a high⁃precision 3D geological model, it has changed the interpretation of traditional horizons to avoid the phenomenon of sandbodies in adjacent layers through smooth interpolation, which ensures that the “stable area” near the fault is as accurate as the one far from the fault zone. The position of the “complex area” of the fault is determined by the coherent properties of the mudstone layer, the range of the “complex area” of the fault is identified by the amplitude attribute, and the inverse attribute interpolation method is used to realize the calculation of the inversion attribute value of the fault crushing zone, thereby reducing the effect of faults on seismic inversion accuracy. Using this method, the accuracy of prediction of sand thickness near the fault was increased from 80.4% to 91.3%, an increase of 10.9%. This method improves the description accuracy of the sand body near the fault and provides a new idea and method for the fine description of sand bodies in the fault attachment.